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Abstract 


Large isotopic anomalies are observed in some abundai t elements in the 
space medium, stellar atmospheres, comets and carbonaceous meteorites. In 
some cases these are due to chemical fractionation, in others it is difficult 
to distinguish between chemical and nucleosynthetic effects. Our research aims 
at experimental and theoretical evaluation of chemical fractionation processes; 
now with emphasis on selective single isotope fractionation in polyisotopic 
systems, particularly in oxygen. The related temperature parameters of % 
meteoritic condensates and of their source medium are investigated by a thermo- 
metric method that is independent of assumptions regarding temperatures and ^ 
pressures in the solar nebula. 

The crucial non-linear chemical fractionation of has been demonstrated 
experimentally in the last year's work under this grant. The effect was 
achieved in condensed CO 2 formed from CO with 12 C 16 0 selectively excited by 
H Lya. The effect was verified by mass spectrometric measurements carried out 
by these isotopic measurements by Prof. R.N. Clayton. 

The meteorite paleotemperature estimates which form a related part of our 
research have been advanced from defining only thermal exposure to evaluating 
time and temperature independently. Grain temperatures at condensation of 
refractory inclusion materials are indicated to be less than 900°K in agree- 
ment with radiation temperature considerations and observations in circum- 
stellar dust shells. 
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1» Scientific objectives 

The research reported here represents a continuing effort to clarify by 
theoretical and experimental investigations what mechanisms are responsible 
for the large chemical isotope fractionation observed by microwave spectrome- 
try in dark interstellar clouds and by laboratory measurements of primordial 
meteorite material >. 

The proposed relationship between the effects seen in the space medium 
and those observed in meteorites is most applicable to the widely varying 
isotopic composition of hydrogen, carbon, nitrogen, oxygen and magnesium and 
possibly some of the effects observed in occluded noble gases, and in odd iso- 
topes of heavy elements such as barium in carbonaceous meteorites. Polyiso- 
topic elements such as oxygen, magnesium and sulfur are of particular inter- 
est in this respect since single isotope enhancement as observed in the labora- 
tory (e.g. , Liuti et al. 1966; Vikis 1978; Turro and Kraeutler 1978) and pre- 
dicted to be effective in the space medium (Haberkcrn et al. 197 6; Arrhenius 
et al. 1979b) can generally be verified only if at least two other isotopes exist 
as reference points. Such enhancement of single isotopes or isotope pairs has 
been shown or predicted to be caused by various mechanisms of selective excita- 
tion and/or transition into reaction complexes of predissociating states (see 
review in Arrhenius et al. 1979b), Selective resonance excitation by spectral 
line sources has been observed astronomically for several molecular isotopic 
species in the space medium (Gahm et al. 1977). Chemical effects of this kind, 
if-found to be the likely cause also of nonlinear isotopic anomalies in meteorites, 
would give useful indications about the physical state of the source medium from 
which primordial solids formed in our solar system. 

An alternative explanation of effects found in meteorites is that not only 

2 16 

small anomalies in heavy elements, but also large effects in, e.g. , H, O, 

1 5 26 

N and Mg, could be of nucleo synthetic rather than chemical origin. Obser- 
vational proof one way or another is generally not possible as far as meteorites 
are concerned, since the effects recorded in them are caused by events in the 
past, and since confirming nucleosynthetic experiments can in most cases not 
be carried out. In contrast, observations of currently active processes in the 


/ 
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interstellar cloud medium provide direct evidence of large scale chemical 
isotope fractionatior in carbon, hydrogen and nitrogen, and with more refined 
techniques the potential exists for discerning the details of selective fractiona- 
tion of individual isotopes of oxygen and magnesium (Winnewisser 1979). Selec- 

l6 17 18 

tive enhancement of O relative to * O in CO has been observed in some 
regions of the interstellar medium and is considered as a puzzling feature if 
interpreted in terms of nucleosynthetic origin (Penzias 1980). 

The chemical isotope effects result from the large thermodynamic 
disequilibria which arc characteristic of the space medium. Since the 
i6otopically "anomalous" meteorite minerals provide a fossil record of the 
conditions in the space medium in the formative era, reliable estimates of 
their formation temperatures are desirable as part of the information needed 
to understand the origin of the isotopic anomalies. 

Speculation about the thermal evolution of the refractory oxide and sili- 
cate minerals containing the oxygen anomalies, including "condensation tem- 
peratures" has been based on observed phase relationships in meteorites 
coupled with hypotheses about the thermal state of the source cloud from which 
the solids formed. Si- ce the outcome of such considerations is critically 
dependent on assumptions concerning temperature and pressure, and since these 
are a priori unknown, the conclusions remain hypothetical. In those cases 
where they are based on equilibrium considerations the conclusions also appear 
unrealistic by an order of magnitude or more when applied to the region of the 
condensation of the most refractory solids (Arrhenius 1971; De and Arrhenius 
1977; Arrhenius and McCrumb 1978). A realistic evaluation of the thermal 
record should preferably be based on measurements of thermally sensitive 
parameters that are' pressure independent. 

Interdiffusion between coexisting solids fulfills this requirement and is 
particularly sensitive in the temperature region down to about one half of the 
melting temperature. Simple diffusion couples such as pure platinum spherules 
overgrown with nickel iron, discovered and measured by El Goresy et al. (1978) 
.were in earlier work under this grant selected for evaluation of thermal exposure 
during and after formation of the nickel-iron overgrowth on the platinum spheres, 
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and of the melilite crystals in which the metal couples are encased. The 
initial choice of pure platinum as an end member avoided the larger uncertain- 
ties about the course of diffusion in the more complex platinum metal alloy 
couples also found in carbonaceous meteorites. 
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2. Theoretical investigations 

Three possible mechanisms 

As pointed out above and dis ussed in detail in the publications result- 
ing from this work, selective isotope effects can arise at molecular excitation 
(e. g. , Liuti et al. 1966), predissociation (Murrell and Taylor 1969: Schaefer 
and Miller 1971: Arrhenius ct al. 1979) or vibrational relaxation (Basov et al. 
1974; Arrhenius 1976). 

Fractionation by transition into repulsive state at predissociation 
Our earlier theoretical work (Arrhenius et al. 1979a, b) specifically exam- 
ined the isotopic selectivity of Franck- Condon type in predissociation. Numeri- 
cal approximation indicated that while predissociation in the classically for- 
bidden region is unlikely to lead to large isotope effects, transition into a repul- 
sive state in the classically allowed region could indeed cause such major 
effects (Fig. 2. 1. 1). Exact calculation proved to be a demanding computational 
task. We decided to postpone further efforts in this direction until experiments 
aiming at isolating this effect from the two others mentioned, have specifically 
indicated it to be quantitatively significant. 

Single isotope selection by resonance excitation 
Previous work 

Much of this year's theoretical work, accompanied by experiments, has 
been devoted to evaluation of this mechanism, both in the laboratory and in the 
space medium (Arrhenius et al. 1980; Corrigan et al. 1980). The fundamental 
theory for the underlying phenomena is relatively straightforward and the achieve- 
ment of highly selective isotope fractionation through this mechanism using laser 
line sources for excitation is well known. The efficiency of UV emission line 
spectra in resonance excitation of CO was first demonstrated by Harteck and his 

associates (Liuti et al. 1966; Dunn et al. 197 3) and later by Vikis (1977,1978) as 

1 3 

discussed above, who achiever, ten and sixfold concentration, respectively, of C 

in product CO and in the central carbon atom of C-O-. A sevenfold concentra- 
18 * 3 c 

tion of O was found in the same photochemical decomposition of CO to CO^ 
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Fig. 2.1.1. Magnitude of overlap integrals (inserts) between the eigen- 
functions for the bound (initial) and repulsive (end) states 
determines the predissociaticn probability. Two isotopic 
states, Ej and E2 of the same vibration are schematically 
indicated, one (above) with large and one (below) with small 
overlap integrals. 
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Ma tchi ng of emiss ion a nd abso rption line a in the laboratory 

and in the s pa ce medium 

The close line matching needed for resonance excitation might 
intuitively seem to make this an exceptional phenomenon achievable only by 
tunable sources or by combining emitters and absorbers selected among a 
large number of species, most of them unlikely to be of importance in nature. 
What makes the phenomenon common in astrophysical settings, however, 

(Gahm et al. 1977) is the band structure of the spectra of abundantly occur- 
ring molecules such as CO or O^; with each electronic -vibrational level is 
associated a large number of rotational levels which may be populated under 
excitation conditions characteristic of the space medium (Fig. 2. 1. 2 ; adapted 
from Tilford and Simmons 197 2 and from Purcell and Tousey I960). This 
together with the occurrence of several strong emission lines in stellar UV 
spectra makes appropriate coincidences sufficiently common for isotopic 
selection to potentially occur in a number of quantitatively important molecu- 
lar reactions in the space medium. Past year's finding of such a coinci- 
1 2 16 

dence between a C O absorption band and the hydrogen Lyman a emission 
line is a case in point (see Fig. 2. ' . 2). The same effect that we invoke here 
as an excitation source for isotopically selective reactions is being used as a 
tool in the detection of specific molecules in the space medium (Gahm 1977). 
The technique is used for abundant molecules that surround astrophysical 
objects or regions that have among others strong atomic emission in the 
Lyman and Balmer lines of hydrogen. The excitation of CO by Lyman a is 
considered to be undesirable for astrophysical detection purposes because 
so many rotational lines are excited in CO (Fig. 2. 1. 2). The more rota- 
tional lines excited the more diffuse the fluorescence in terms of spreading 
the input energy as output over many more lines. This astronomically 
detrimental circumstance is advantageous to our work because we are inter- 
ested in efficiency of the molecular reactions which take place before 
fluorescence can occur. 

Of the astrophysical settings that Gahm et al. consider, an interest- 
ing one in terms of explaining oxygen effects in meteorites is the circumstcllar 
region around emission line stars of spectral classes O through G, and the 
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Hydrogen Lyman a Emission from Sun 


Hydrogen Lyman a Emission 
from typical laboratory source 
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Fig. 2. 1. 2, (Below) Overlap of laboratory Lyman a emission line with rota* 
tional lines of CO in the (14*0) vibrational band of the A-X 
electronic transition. 

The figure illustrates that emission overlaps absorption lines 
of ^-^C^O-isotopic species only, leading to selective excitation 
of this species. 

(Above) Expansion of scale of 1 C O portion of (14-0) band with 
schematic rotational absorption lines and with Lyman a emission 
profiles from laboratory sources and from the Sun. Emission pro- 
file width is to scale relative to location of and l2 C^O- 

bandheads. Absorption lines of other isotopic CO species are off 
the figure to the red. 

The broadening of the stellar emission profile (with overlap 
over about 40 rotational l^C^O lines) leads to an order of mag- 
nitude enhancement in astrophysical settings over the effect 
achieved in the laboratory where the narrow profile lamp source 
can be made to effectively overlap only two such rotational lines. 
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T-Tauri stars. The Sun is now a G2v star and in its early evolution it is 
believed to have passed through a T-Tauri stage. Several other situations 
where dense proto stellar clouds are irradiated by strong line sources have 
also to be investigated. 

Consideration of such selective absorption also forms the baMs for 
our studies of resonance excitation of isotopic lines in the Schumann- Runge 
bands of and in the fourth positive band of CO, both important absorbers 
in space and that of also in the Earth's atmosphere. This work led to the 
discovery (McCrumb 1979) of a secondary source of ozone in the Earth's 
upper atmosphere (Cicerone and McCrumb 1980). 
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3. Experimental studies of selective isotope fractionation 
Background 

Laboratory experiments guided by the theoretical considerations above 

provide a means for clarifying the mechanisms in space that lead to the large 

isotope fractionation observed to occur there. Experiments carried out under 

the present grant have demonstrated large isotope effects in reactions of 

excited molecules and c Secular ions in the systems C-0 ar.d N-O. Fraction- 

15 13 

ation up to 42 per cent was observed in N (Fig. 2. 2. 1), 5. 6 per cent in C 

18 

and 6. 2 per cent in O (Fig. 2. 2. 2), extended to 10 per cent in earlier work, 

and ~ 50 per cent in last year's UV experiments (Table 2. 2). 

The oxygen system is of particular interest in view of the substantial 

enhancement of found and explored in meteorites by R. N. Clayton and 

his collaborators. Much of our efforts during the last few years for this 

reason have been devoted to studying the fractionation of the tl\ree oxygen 

* 

isotopes in reactions of excited molecules. The excited system (C-O) was 
used for these model studies because of our earlier findings of large isotope 
effects, the importance of the system in protostellar cloud chemistry, ex- 
perimental convenience and the experimental demonstration of large selec- 

18 

tive fractionation of O (sevenfold) in product CO^ and C^O^by resonance 

VUV excitation (Liuti et al. 1966; Vikis 1978; further references in Arrhenius 

13 

et al. 1980). Concurrent fractionation of C was also monitored in our 
investigations. 

Selective enhancement of in product CO , from H Lyq 
irradiation of CO 

Much experimentation was during the last year devoted to increasing 
the fractionated product yield from the H Lya stimulated CO reaction. It was 
assumed (as later supported by observation) that as a first approximation the 
selective fractionation of * O would be proportional to the linear fractionation 
in the product CO^. The latter generally exceeded a few per cent and could 
be routinely monitored with the SIO Micromass spectrometer. As a result 
of various technical improvements the yield could be raised sufficiently to 
permit the production of samples large enough (10-20 |imoles) to be used by 
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2. 1. Nitrogen isotope fractionation at formation of nitrogen 
oxides from a cold plasma (from Arrhenius et al. 
1978). 



2. 2. Fractionation of and l *C in excited state reactions 
at different kinetic temperature, pressure and flow 
rate (from Arrhenius et al. 1979b). 
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Prof. R. N. Clayton for fluorlnatlon and Isotopic analysis; due \o his heavy 
commitments only two sasples have been analysed at the present tine. Both of 
these show anomalies of the type found in refractory inclusions and of a mag- 
nitude that makes Prof. Clayton convinced that they are real. These results 
must, however, for reasons discussed earlier, be Interpreted with caution 
until the isotopic composition of the source CO has been verified. This is 
one of the most Important undertakings for the completeness of demonstration 
of 16 0 single isotope fractionation by resonance excitation. Similar frac- 
tionation by isotopic filtering has already been demonstrated by Bally and 
Langer and has in last year's work under this grant been applied to the Earth's 
atmosphere. The results indicace that this effect provides a second source 
of ozone below the present ozone maximum (Cicerone and McCrumb, 1980). 
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Figure 2. 2. 3. Oxygen isotope fractionation data from R. N. Clayton 
(see also, table above). 
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To minimize the effects of all types of backgrou i gases (hydro- 
carbons, H O, CO^) all seals are metal (couplings and valves) andbakeable 
to 300 a C. The stainless steel reaction vessel (12 cm dia. X 35 cm long; 
Figs. 2. 2. 4 and 2. 2. 5) is coupled to the turbo pump by a 3. 4 cm diameter 

vacuum V ae incorporating a bakeable metal sealed valve. This line is used 

-8 

between experiments to evacuate the reaction vessel to 5 X 10 torr during 
bake out. Bake out times are typically 12 hours. 

The reaction vessel lid, incorporating the UV lamp body and a 
25 X 1/2 mm MgF^ window, is sealed to the reaction essel by a 6. 25 cm 
copper gasket. 

The lamp body consistr; of a 9 X 12 mm quartz tube 12. 5 cm long 

with one end fused to a stainless steel tube of the same diameter. The gas 

flowing through the lamp is excited by an air-cooled Evenson cavity powered 

by a 2. 5 GHz magnetron. The lamp gas flow system has a gas supply bottle, 

pressure regulator, needle flow valve, pressure gauge (1000 torr), and LN 

trapped vacuum waste line (Fig. 2. 2. 3). The lamp can be operated in a 
% 

static mode or with continuously flowing source gas. For optimum emission 
of H Lya a gas mixture of 2% H , in He was used. The 2. 5 GHz lamp is 
typically operated at 40 watts with a gas pressure of 1 torr. 

The sample handling system consists of a bakeable stainless steel 
trap, two digital pressure gauges (0-1 and 0-1000 torr) and gas sample 
extraction port. This section is connected to the reaction vessel via a four- 
way cross, an integral part of the 3. 4 cm diameter pump line. 

The quadrupole mass spectrometer used for on-line mass spectrome- 
try is used to monitor impurity levels in the VUV source and in reactant 
gases. This spectrometer has also proven most useful in monitoring gas 
reactions as a function of pressure and temperature, but with limitations 
that will be eliminated with the readout system requested. For determina- 
tion of variations smaller than about 20 per cent, the line to background ratio 
is, however, limiting in the present system; for all such measurements we 
rely on magnetic sector mass spectrometers. 
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Figure 2. 2 . 4. Stainless steel reaction vessel, feed lines for reactant and VUV 
lamp gases, gas sample handling, and mass spectrometer sys- 
tems for VUV resonance excitation experiments. The entire 
system is stainless steel, bakeablc to 300 *C with an ultimate 
vacuum pressure of 1 X 10"® mm of Hg. 









A. Reaction vessel 

B. H Lyman a lamp 

(j i Air-cooled Evenson cavity 
D. Turbomole cular pump 


Figur. 2.2.5. App.r.tu. for ..Uctiv. r«.o«nc« .experiment, 
schematically shown in Fig. 2. 2. 4. 
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4. Thermal exposure of grains with enhanced O content 
Background 

In the space medium, as observed today, there are typically large 
differences between rotational, vibrational, translational, neutral, ion and 
electron temperatures. High internal molecular temperatures coupled with 
low kinetic temperatures give rise to the large isotope effects maintained 
in steady state in chemical reaction cycles. Isotopically selective resonance 
excitation is commonly observed by fluorescence; we have proposed that 
this excitation results in isotopically enhanced reaction products including 
refractory oxides in meteorites. 

The way in which gas and plasma temperatures couple to solid 
particles determines the thermal state and growth conditions in condensing 
grains. We wish to reconstruct these conditions for meteoritic grains using 
isotopic anomalies as a means for understanding the state of the source 
medium from which they formed, and as built in solid state thermometers 
for evaluating the temperature of the grains themselves. 

/. fundamental requirement for condensation and growth to proceed 
is that energy is removed frcm the system, primarily by radiative cooling 
of the grains since solids are effective IR emitters. Radiative cooling can 
only take place if the medium surrounding the condensate is optically thin. 
The rate of condensation in contracting dense regions of protostellar clouds 
should consequently be inversely proportional to the optical depth of the 
locale. This has important consequences for the temperature differential 
between the condensing grains and the source medium; the solid particles 
characteristically maintain steady state temperatures an order of magnitude 
below the kinetic temperature of the medium. Accordingly, in the — 7 000° 
"cool" plasma in nova shells dust is observed to emit at ~ 1000°K (Geisel 
et al. 1970; Clayton and Wickramasinghe 1976); dust particles embedded in 
60, 000* plasma in the Jovian magnetosphere assume temperatures below 
200* (cf. , Fig. 2. 3. 1). The radiation balar.ce determining the temperature 
differential between grains and source medium at condensation has been 
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obKrvcd phenomena 



Figure 2. 3. I. The astrophysical P-T setting for chemical reactions, 
isotope fractionation and condensation in interstellar 
and protostellar clouds with embedded or surrounding 
UV line radiation sources (from Arrhenius and Raub, 1978). 


calculated by Lehnert (1970),. Arrhenius (1971) and, as a part of work und^r 
this grant, by De and Arrhenius (1979). 

Under these circumstances, it is impossible to determine "condensa- 
tion temperatures" even within an order of magnitude from assumed thermo- 
dynamic equilibria, and without distinguishing the relative contributions from 
isobaric cooling and isothermal compression. The only way to reconstruction 
of actual grain temperatures would seem to be by measurement of such 
phenomena in the grains which are sensitive to their thermal exposure, i. e. , 
the length of time which they have spent at specific temperatures. 

The intermetallic diffusion studies carried out under the present grant 
give stringent upper limits for the temperature, and less precisely for the 
duration of thermal exposure of the refractory minerals which are the carriers 
of some of the observed isotopic anomalies, notably in oxygen. 
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The picture that emerges from these studies is an initial condensation 
of the most refractory oxides and silicates at a grain temperature which could 
not have exceeded ~ 1000*K and which could have been much lower (Arrhenius 
and Raub 1978; Arrhenius, Raub and Schimmel 1980). The excitation tempera- 
ture (in contrast to the kinetic temperature) of the source medium, would have 
been sufficiently high to achieve the observed isotope fractionations in nitrogen, 
carbon and hydrogen (not necessarily in the same region and at the same time), 
and at least the mass dependent effects in oxygen and magnesium in carbon- 
aceous meteorites. 

Earlier measurements 

Using diffusion constants and thermal activation energies for Pt-Fe 
and Pt-Ni in earlier literature, we were able to demonstrate (Arrhenius and 
Raub 1978) that the melilite and spinel, in which the Pt-(Ni, Fe) metal couples 
occur abundantly, could only for a small fraction of a second sustain the tem- 
. perature of 1625 K inferred from chemical equilibration considerations and 
an arbitrary pressure assumption to be the "condensation temperature" at 
which the mineral assemblage formed (Grossman 197 2). A grain temperature 
of 900 K could be sustained for a day and 700 K for a year; growth times of 
this order probably represent the minimum necessary for the formation of 
the refractory grains investigated. 

The phase relationships, including ordering, in the systems Fe-Pt and 
Ni-Pt are simple and well-known, and the interpolation of phase structure and 
diffusion kinetics of Pt-(Feg ^Ni^ ^), the composition observed in carbon- 
aceous meteorites, appears reliable at the level of accuracy needed for the 
discussion above. Larimer (1979) proposed that unknown phases may 
occur in the ternary system, which would extensively slow down inter- 
metallic diffusion so as to invalidate the conclusions from the experimental 
data quoted above. The only phase formation (at low temperature) or precursor 
phenomenon to be expected in this system is, however, due to ordering at stoichio- 
metric compositions. As is well known and also shown for this system in our work 
(see below) they result in acceleration of the diffusion rates, not in slowdown. 

An experimental study was undertaken during the last grant year 

of the diffusion of Pt into Fe. .Ni n . at a series of temperatures, and the 

0. 4 II. I* 
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diffusion coefficient determined ss a function of composition in the diffusion 

path (Arrhenius. Raub and Schimrnel 1980). The diffusion coefficients thus 

-6 2 

measured in this ternary system at 1100 K range from 4 . 6 * 10 cm /day 
at 5% Pt to 7 • 10‘ 6 at 10% Pt, in concord with earlier published data for the 
binary systems at the same temperature. For these Kubaschewski and 
Ebert 0944) found 1.0 . 10~ 5 for Pt Q g 5 Ni Q 15 »nd Kucere and Million (1975) 

4. 0 . 10 ^ for PIq ^Fe^ 4 and 8* I * 10*^ cm^/day for Pt Q ^Fe^ g* 

The effect of diffusion observed at this rate for one day at 1100 K was 
to lower the platinum content in the original 5 |lm thick pure platinum layer 
of the diffusion couple to 23 per cent, thus greatly exceeding the diffusion 
alteration of the platinum spherules in the meteoritic melilite- spinel, which 
is indistinguishable from 100% Pt. 

These measurements thus confirm the conclusion that the melilite- 
spinel-Pt-FeNi assemblages cannot have formed at grain temperatures 
exceeding about 1000 K or have been subsequently exposed to temperatures 
of the solid in this range for more than a few minutes to an hour. Formation 
of the Pt-(Fe,Ni) couples by exsolution is obviouslv precluded by the unlimited 
mutual solid and liquid solubility of these metal; . 

The physical basis for the relatively low temperature of the condensing 
refractory solids indicated by these measurements is well understood from 
the astronomical observations and related theoretical considerations discussed 
above. Condensation can be brought about by the observed densification of 
the dilute cold interstellar cloud medium (gravitational or hydromagnetic con- 
traction) or, more hypothetically, by cooling of a hot, dense cloud, or by a 
combination of these. Each case or combination would result in practically the 
same condensation sequence; subtle indicative differences are interesting to 
exploit. 

If the temperature of condensing refractory grains were sufficiently 
low relative to the activation energy for crystallization (100-200 K) a refrac- 
tory glass would result. The high latent heat of crystallization of the stable, 
crystalline compounds, released by collision or radiation heating, is capable 
of momentarily bringing grain temperatures up to the order of 1000 K, result- 
ing in the observed "igneous" texture (Arrhenius 1976). 


21 


ORIGINAL PAGE 13 
OF POOR QUALITY 

Oxide diffusion barriers, which have also been invoked as retardants 
arc thermodynamically excluded in this system at temperatures as high as 1600 
1700 K and at the reducing gas compositions generally assumed. Measurement 
of diffusion of metals through iron oxide films demonstrated them to be in> 
effective as diffusion barriers at the thickness that would have to be involved 
(less than visible by scanning electron microscopy) (Cabrera 1980). Further- 
more. special and unlikely conditions are required for an oxide layer to form 
on platinum spherules, followed by growth of metallic nickel -iron. 

Information from diffusion experiments during the last grant year 

The investigations summarized above confirm that diffusion measure- 
ments in platinum -nickel iron couples overgrown with refractory silicates in 
carbonaceous meteorites provide observational evidence for the upper limits 
of temperatures for which these refractory minerals have been exposed during 
their lifetime, including their formation and metamorphism. Considering 
the still more detailed information potentially retrievable from these and 
different types of diffusion couples in meteorites, a series of new experi- 
ments were undertaken. They aim at measuring the difference in diffusion 
rate between iron and nickel in platinum and utilizing this difference as an 
independent temperature indicator. With measurements of both total and 
differential diffusion it should be possible to separate the time and tempera- 
ture parameters in the fossil diffusion record. 

Theoretical considerations 

Diffusion theory is based on two assumptions. First that the rate of 
transfer of matter, the flux, across a plane perpendicular to a chemical gradi- 
ent is proportional to that gradient. This is known as Fick's first law and is 
written F - M • Vu. Second, it is assumed that matter is conserved over time, 
t, for a closed system, i, e, 7F + 7C = 0. 

For the case of diffusion in a binary alloy along a unidirectional composi 
tion gradient the flux, 

J = -D • dC/dx (1) 

2 

where J is in units of g/cm / sec, and'D, the diffusion coefficient, is in units 

’ • 

of cm^/scc. If we choose a fixed reference plane, thus fixing the diffusion of 
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both components to one coordinate axis, then B, the interdiffusion coefficient, 
measures both the rate of movement of the individual atoms and the bulk flow 
of the lattice (Kirkendall effect), caused by differences in the rates of diffusion 
of the components. 

Assuming a closed system we derive a relationship which can be dealt 
with experimentally, namely 

dC/dt = -d/dx • (-D • dC/dx) (2) 

For most systems the Interdiffusion coefficient is a function of concen- 
tration. Matano (19 ) developed a solution assuming that certain conditions 

are maintained during the experiment: (1) for t = 0, C=0 at x<0, and 
C: 1 at x> 0, (2) for all t, dC/dx], = 0, and (3) xfjt = C, where C is 
constant for a given concentration, x the distance migrated from the 
"Matano interface" and t the time period of migration. By substituting 
X = x/7T into equation (2), rearranging and then integrating, it is found that 

r c 

/ xdC = -(2t)(dC/dx)S 
0 C 

Given a c-x curve for some time temperature combination, we can calculate 

, where c denotes a specific composition. 

It is found that the temperature dependence of B is typical for 

c 

chemical reaction rates in that it follows the expression 

D . (D ) exp(Q /RT) 
c o c c 

The pre- exponential factor (D ) and the activation energy Q , which is re- 

o c c 

lated to the bond enthalpy, can be determined by graphing 2a vs 1 /T. 

Diffusion in an n- component system is described by an extension of 
Fick's first and second laws (Ktrkaldy 1958). The flux of a component i, 
relative to a fixed reference plane is given by 

n - 1 

8 ‘ 2 8" dC./dx , i * 1, 2, . . . , n-1 
1 j = 1 4 j 
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where D?, is an interdiffusion coefficient with n as the dependent variable, 
and dCj/dx is the gradient of the component j, there being n-1 independent 
concentrations. By Fick's second law 


n-1 

dC /dt = d/dx • X S” dC /dx 
1 j*l 3 J 

The Matano evaluation can be extended to the multi -component system 
(Kirkaldy, 1958) and the solution for By is found to be 




xdC. 

l 


n-1 




dCydx 


For the ternary system of Pt-Ni-Fe this reduces to 


/ xdC [®N iNi dC Ni ,dx + ®NiFe 

C Ni 

o 

f Fe X dC = - 2t f B ^ dC /dx + D 1 * 4 
j xqv, « |_FeNi Nr FeFe 

C Fe 


d Cre /dx] 

dC Fe / dx ] 


(Note that any element can be taken as the dependent variable. ) Since each 
of these equations contains two unknowns, it is necessary to design the experi- 
ment so that, for a given temperature, two separate diffusion paths will inter- 
cept simultaneously for all three components. In most studies the initial 
diffusion couples cover a wide range of compositions. 

Assuming that these ternary diffusion coefficients also exhibit 
normal reaction rate behavior, then for a given diffusion path evaluation for a 
specific concentration requires the simultaneous solution of four equations. 
Since in virtually every system the species are diffusing at different rates 
it is expected that the solution would yield a unique temperature, a result of 
primary importance here. 
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Experimental and analytical procedure 

The initial semi-infinite experimental media used last year simulated 
typical natural assemblages of pure platinum spherules in Ni^Fe (Arrhenius 
and Raub 19? 8); the diffusion couple consisted of a 5 pm thick layer of pure Pt 
sputtered onto a sheet of Ni^Fe, 1 mm thick. 

Thic year's experiments employed an infinite source for both the Pt 
and the Ni-Fe which yields more extensive information about diffusion proc- 
esses in this system. The Pt was prepared from Pt foil (with impurities of 
less than 1%). which was arc melted under argon atmosphere into a bead. 

This was pressed into a disc, then the surfaces were ground flat and polished 
with 1 pm diamond grit. The discs were annealed at 1100" for 5 days to 
develop the crystals and to remove imperfections in the lattice. Then the 
surface of each disc was lightly repolished with 1 Mm diamond grit to remove 
any contamination. 

The Ni^Fe from the previous experiments was utilized, the plates 
were poli-hed in the same manner as the Pt discs and annealed at 900* for 
5 days. Again, the surfaces were polished with 1 pm grit. The diffusion 
couple 8 were fixed in a small vise of stainless steel and were sealed in quartz 
tubing under argon atmosphere. The couples were welded for 1 hour at 900", 
quenched and then sectioned. The diffusion experiments were then carried 
out at temperatures ranging from 900"C to 1300*C, with an error of ±5*. 
They were then mounted in epoxy and prepared for quantitative analysis. 

The concentration measurements were done on a Cameca electron 
microprobe, and were carried out in two stages. Each couple was scanned 
along numerous traverses while recording the analog, and in the quantitative 
measurements the digital output for each of the elements. All results were 
found to be reproducible within ± 1%. 

Results 

Eleven samples were annealed at 117 3*K (68.5 and 183 hours), 1223"K 
(74. 5 and 121 hours), 1273" (39. 5 hours), 1373", 1473" and 1573"K (3. 99 and 
5. 53 hours). Three concentration profiles are presented in Figs. 2. 3. 2a - c. 
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Caption for Figure 2. 3. 2 -- Next page. 


Diffusion profiles for Fe and Ni in Pt, and Ft in Nl^Fe at three 
different temperatures (three graphs) and after different lengths of 
diffusion time (indicated in graphs; 16. 2 hours at 137 3*K). 

Two interesting features can be seen. One is the beginning of order- 
ing at 1223*K, manifested by acceleration of the diffusion in the composi- 
tion region around Pt Q< jFsq^ 25 Ni 0. 25 (Ptl**eNl]). 

The other is the rapid increase in the diffusion rate of iron rela- 
tive to that of nickel, beginning between 1223 and 137 3 *K. This phe- 
nomenon provides a temperature indicator independent of diffusion time. 

The degree of interdiffusion of nickel in platinum in graph c), for 
example, can be achieved during 16. 2 hours at 137 3°K as shown, but 
also at an infinite number of other combinat’ons of longer times with 
lower temperatures and vice versa; at higher temperature the times 
rapidly -become unrealistically short. A distinguishing feature is, however, 
that Fe is seen to have raced ahead of Ni into the Pt, pinpointing the tempera- 
ture as being above 1280*K. On the basis of this and the total diffusion, 
the time is fixed at around 1 day or less. 

In the meteorite spherules investigated no perceptible diffusion of either 
iron or nickel has taken place. 
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Figure 2, 3, 2 




33 


ORIGINAL PAGE IS 
OF POOR QUALITY 

The position of the initial boundary between the members of the diffusion 
couple was in each sample indicated by a plane which contained voids 
developed from the separation of minute argon pockets during annealing. 

All analyses were made in regions without voids. 

Any ambiguity about the position of the initial phase boundary in the 
diffusion couples does not affect calculation of the diffusion coefficients, 
because movement is taken relative to a hypothetical plane, the "Matar.o inter- 
face, 11 where the diffusion is equivalent in both directions. Table 2. 3. 1 con- 
tains the interdiffusion coefficients for the pseudo-binary system Pt-M, where 
M = Fe, Ni. The assumption that Fe and Ni behave like species is reasonable 
on the basis of similar phases in the Ni-Pt and Fe-Pt phase diagrams. Also, 
Fe and Ni form a continuous solid solution and there is only one ordered 
phase, Ni ^Fe, which has a low thermal stability. Our continued work aims 
at making the additional still more elaborate measurements that will be 
needed to evaluate the ternary diffusion coefficients. In the present work 
we are interested in observing the system over a broad range of temperatures. 

* Discussion 

It. is apparent from the concentration profile that, as expected, ordered 
phases occur in tnis system; corresponding irregularities in the curves are 
particularly distinct for 117 3*, 1223* and 127 3*. When translated into diffu- 
se, rates (Fig. 2. 3. 3), this ordering is represented by maxima and minima 
i : o>xnd the stoichiometric compositions MPt^ and MPt. MPt is the most 
stable and is the only phase to affect rates (by accelerating them) at 1000*. 
These observe ’°ns correlate well with the known ordered phases in the Ni-Pt 
and Fe-Pt sys ms. They all show the strongest ordered phase at 50% Pt. 

In the binary systems ordering is observed at 25% Pt while in the ternary 
system ordering is more pronounced at 75% Pt, where also the phase FePt^ 
has been reported. 

Indications that ordering can affect diffusion at temperatures well above 
the order -disorder transition has been presented before, Kucera and Million 
(1975) used the phenomenon of short range ordering to explain anomalies in the 
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Table 2. 3. 1 

Interdiffusion coefficients for the pseudo-binary system Pt-(Ni, Fe) 
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* 


10 *0 50 *0 Ho 



Figure 2. 3. 3. Diffusion constant as function of composition at various tem- 
peratures, showing the increased rate of diffusion around 
integral stoichiometries (25, 50, 7 5 at. % Pt) and at low 
temperatures where ordered compounds begin to form. 

In contrast to what has been implied by Larimer (1979). the 
formation of ordering compounds in this system thus does not 
impede but accelerates diffusion and makes the present tech- 
nique even more sensitive in these composition -temperature 
. regions. 
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diffusion of Ni into NiFe. The rate of diffusion peaked at the composition 
Ni^Fe for diffusion temperatures of 814*C, 924®C and 1076 *C. The actual 
phase transition occurs only below 503 *C. They also interpreted data for 
diffusion of Co into NiCo in a similar way. The rate peaked at NiCo^ and 
Ni^Co for temperatures below 1300®C. The ordering transition for Ni^Co 
is reported to occur at and below 1033 *C. 

In all the concentration profiles in this study (e.g. those shown in 
Fig. 2. 3. 2), it can be seen that Ni^Fe is a particularly stable composition. 
This is most apparent for 117 3° (68. 5 hours)(Fig. 2. 3. 3). Yet the actual 
formation of this phase occurs only at and below 77 3®K. An important point 
is that the ordering phenomena result in acceleration of the diffusion rates, 
making our conclusions about thermal exposure conservative (Arrhenius and 
Raub 1978; Arrhenius, Schimmel and Raub 1980). 

Also evident from tne concentration profiles is the interesting and 
useful fact that the relative rates of diffusion of Ni and Fe change with tem- 
perature. For 127 3° and below Ni penetrates more deeply than Fe; for 
137 3* and above the reverse is true. It is interesting to note that tracer 
diffusion data by Kucera (1975) for Pt into Pt-Fe shows a large jump in the 
rate of diffusion for temperatures above 127 3°. 

Figure 2. 3. 4 contains graphs of Bn against l/T. If a single rate 

were involved, these points should lie on a straight line with slope Q /R. 

c 

Obviously, this is not the case and we choose to interpret the data as repre- 
senting a break in diffusion rates between 127 3° and 137 3®K. Irregularities 
below 127 3° are due to the short range ordering preceding the formation of 
ordered phases. Above 127 3® , the rapidly increasing disorder leads to a 
change in bond energy which affects the diffusion rate of iron more strongly 
than that of Ni. 

A major conclusion from these results is that in multicomponent diffu- 
sion systems of the kind occurring in meteorites, the temperature dependence 
of the diffusion rates of individual components such as Fe and Ni permits us, 
in principle, to determine the fossil temperature of diffusion as a parameter 
separate from the total extent of diffusion, which measures total thermal 
exposure (time X exp. temperature). 
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Caption for Figure 2. 3. 4 -- Next page. 


Plot of the (negative) logarithm of the diffusion coefficient against 
inverse temperature for increasing platinum concentrations in the system 
Pt-Fe, Ni (top row 10, 20 and 30% Ft; second row 40, 50, 60; third 70, 

80, 90% Pt). The diagrams show over this entire concentration range an 
increase in the diffusion rate with onset around 1323*K (1050*C). The 
increasing dependence of the diffusion constant on stoichiometry as tem- 
peratures drop into the ordering range is also seen, with anomalously 
high diffusion constants for the incipiently ordered compounds (cf. also 
Fig. 2.3.2). 
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